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ABSTRACT 

We have conducted a 3-month program of simultaneous optical, soft and hard X-ray 
monitoring of the LMXB 4U1636-536/V801 Ara using the SMARTS 1.3m telescope 
and archival RXTE/ASM and Swift/XRT data. 4U1636-536 has been exhibiting a 
large amplitude, quasi-periodic variability since 2002 when its X-ray flux dramatically 
declined by roughly an order of magnitude. We confirmed that the anti- correlation 
betwe en soft (2-12 keV) and hard (> 20 keV) X-rays, first investigated by IShih et al.l 
(2005), is not an isolated event but a fundamental characteristic of this source's vari- 
ability properties. However, the variability itself is neither strictly stable nor changing 
on an even longer characteristic timescale. We also demonstrate that the optical coun- 
terpart varies on the same timescale, and is correlated with the soft, and not the hard, 
X-rays. This clearly shows that X-ray reprocessing in LMXB discs is mainly driven by 
soft X-rays. The X-ray spectra in different epochs of the variability revealed a change 
of spectral characteristics which resemble the state change of black hole X-ray binaries. 
All the evidence suggests that 4U1636-536 is frequently (^monthly) undergoing X-ray 
state transitions, a characteristic feature of X-ray novae with their wide range of lumi- 
nosities associated with outburst events. In its current behavioural mode, this makes 
4U1636-536 an ideal target for investigating the details of state changes in luminous 
X-ray binaries. 
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1 INTRODUCTION 

The low mass X-ray binary (LMXB) 4U1636-536 is a well- 
known type-I X-ray burster, which has been extensively 
observed over the past four decades and is considered 
we ll-established as a persistent X-ray source (see references 
in iLiu et all 120071 ). Its optical c ounterpart (V801 Ara) 
revealed its 3.8h orbital period (|Smale fe Mukai I Il988l ) 
through photometric variations , but the period h as also 
been detected spectroscopically l|Casares et al.|[2006l ). How- 
ever, during 2001, the All Sky Monitor onboard the Rossi 
X-ray Timing Explorer (RXTE/ASM) recorded an obvious, 
gradual decline of its X-ray flux level. Within a year, it 
began to exhibit a substantial quasi -periodic, long-term 
variability of ~46d l|Shih et all 120051 ). More remarkably, 
Shih et al. also reported that this variation was present 
at both soft (2-12 keV, RXTE/ASM) and hard (20-100 
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keV, INTEGRAL/IBIS) X-ray energies, but that these 
were anti- correlated. Such a pronounced energy-dependent 
phenomenon led to the suggestion of an accretion disc 
instability as the most likely explanation, as the result of a 
decrease in the mass transfer rate. 



iBelloni et all l|2007f) provided further evidence regard- 
ing the nature of the long-term variability of 4U1636-536 
based on a ~1.5 yr RXTE/PCA monitoring campaign from 
March 2005 to September 2006. The lightcurves display 
long-term variability of 30-40d which is consistent with 
the temporal and energy changes that occur between two 
states o f atoll sources: the island ( hard) and banana (soft) 
states (|Hasinger &: van der Klis I 1 19891 ). From recently 
published Swift /BAT hard X-ray (2 0-50 keV) light curves 
that began in 2005, iFarrell et al] l|2009l ) have confirmed 
the anti-correlated relationship between the soft and hard 
X-rays in 4U 1636-536, providing further support for the 
suggestion of a state change/accretion instability within 
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the source. 

According to the standard model for LMXBs, 
their optical emission is dominated by X-ray reprocess- 
ing in the accretion disc arou nd the compact object 
( van Paradiis fc McClintock|[T99i ). Thus, the varying X-ray 
flux of 4U1636-536 should also influence the apparent 
brightness of its optical counterpart. With this in mind, 
we began a ~100d optical daily monitoring program which 
was expected to cover at least 2 cycles of the previously 
observed long term X-ray variability in 4U1636-536. 

In this paper, we report the results of our optical 
observations, and compare them with the simultaneous 
X-ray data available from the RXTE/ASM and Swift/BAT 
archives. These X-ray monitoring data also allowed us to 
investigate the long-term variability properties of 4U 1636- 
536 as a function of X-ray energy. An X-ray spectroscopic 
analysis of RXTE/PCA and HEXTE data was carried out 
to investigate the evolution of its spectral properties as a 
function of the long-term variability. We also discuss pos- 
sible physical mechanisms that might drive such phenomena. 



2 OBSERVATIONS 

2.1 SMARTS optical campaign 

We performed a 3 months optical monitoring of 4U1636- 
536 using the NOAO SMARTS 1.3m telescope with the 
ANDICAM on Cerro Tololo, Chile between 01 April and 
30 June 2008. The CCD for the ANDICAM is a Fairchild 
447 2048x2048 CCD with 15-micron pixels. Each night, a 
pair of exposures of 2 minutes each through both B and V 
filters was taken. In total, we obtained 67 nights of good 
data during the campaign. Data reductions, such as bias 
subtraction and flat-field correction, were provided by the 
Yale SMARTS teanfl 

We utilised the aperture photometry package 
IRAF/DA0PH0T to extract photometric information from 
the data. A number of bright comparison stars (listed 
in the NOMAD catalog) on the same CCD frame were 
selected to reconstruct the variability of the source, but 
absolute calibration of the field was not performed in 
this work. From the scatter of the differential magnitudes 
of companion stars, the accuracy of the differential mag- 
nitude measurement of 4U 1636-536 is typically ~0.006 mag. 



2.2 X-ray monitoring data 

Since the 1960s 4U1636-536 has been frequently observed 
by many X-ray missions. The All-Sky Monitor (ASM, 2-12 
keV) onboard the long-lived Rossi X-ray Timing Explorer 
(RXTE) has been providing the longest continuous coverage 
of all luminous X-ray sources, stre tching from 1996 to the 
present day fsee lLevine et al1ll996T ). Additionally, the Hard 
X-ray Transient Monitor program, using the Burst Alert 



Telescope (BAT, ^1 5-50 keV) onboard Sw ift, has been 
operating since 2005 ijBarthelmv et al.ll2005l ). We retrieved 
both datasets up to 09 February 2009 from the HEASARC 
archival and the website of the Swift/BAT Hard X-ray 
Transient Monitoi0, respectively. 



3 ANALYSIS 

Figure [1] displays the three months of SMARTS monitoring 
in the optical B and V bands, together with the simultane- 
ous soft (RXTE/ASM) and hard X-ray (Swift/BAT) data. 
The first feature to note is the presence of a nearly regular 
variation in all wavelengths, which is even more obvious in 
the X-ray data. In figure [2] the Lomb-Scargle periodogram 
reveals a significant, consistent period of ~30 d at all 
wavelengths except the optical V band. We notice that the 
peak profile is rather wide, so we simulated our data with 
a 30d sine wave (plus Poisson noise) using a similar count 
rate as observed for 4U 1636-536. This simulation shows 
a peak with similar width to that observed, suggesting 
that it is a result of the short length of the simultaneous 
dataset (only covering a few cycles). Indeed, the peak 
width becomes narrower as we increase the duration of 
data included (RXTE/ASM and Swift/BAT only). For 
example, within one year surrounding the campaign, the 
source clearly exhibited a significant long-term variability 
of ~ 30d in both energy bands. 



3.1 Longer term variability 

The period present in the 2008 camp aign was clearly sh orter 
than that exhibited in 2004 (~46 d. IShih et ail 120051 1. and 
was not always con sistent with other reports, i.e. 30-40d 
l|Belloni et al.l 120071 ) and ~46d l|Farrell et al.1 120091 ). Until 
recently, such variability studies have been rather patchy, 
and the correlation between soft and hard X-rays was not 
fully established. However, with the availability of both 
RXTE/ASM and Swift/BAT monitoring over extended 
periods of time, a time-depe ndent power sp e ctral analysis, 
similar to that reported by IClarkson et al.l (|2003l ) on the 
superorbital modulation in SMC X-l, was performed in 
order to investigate the long-term variability evolution in 
4U 1636-536. As the timescale of interest is ~50d, a sliding 
data window of 200d was used. In order to provide some 
smoothing of the resulting dynamical power spectrum, each 
subsequent window moved forward by 50d. The resulting 
power spectrograms from both datasets are shown in figure 

El 

Our analysis suggests that even the dominant period- 
icity (the darkest region) is transient in nature. The most 
recent (2009) data show that the variation has evolved to a 
shorter timescale (~30d), but our analysis did not reveal 
any ultra long-term periodic tr end, as was seen in e.g. 
SMC X-l (|Clarkson et al.ll2003h . Near the three months 



1 http://www.astro.yale.edu/smarts/ANDICAM/ 
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Figure 1. From top to bottom, Swift/BAT (hard X-ray), RXTE/ASM (soft X-ray) count rates and SMARTS 1.3m telescope optical (B, 
V) lightcurves of the LMXB 4U1636-536. Relativity magnitudes are presented in both optical B and V bands, and the offset is ~17.4. 



of SMARTS observations (HJD ~54600), it appears that 
the ~30 d period was only significant in Swift/BAT data. 
However, as we have shown in figure [2] the ~30 day period 
is also significant in RXTE/ASM data. 



3.2 Correlations 

To investigate the correlation between different wavebands, 
we performed a cross-correlation analysis with the three 
datasets: (i) optical and soft X-rays, (ii) optical to hard 
X-rays, (iii) soft and hard X-rays. Since the variability is 
more significant in optical B band, that was used in the 
analysis (after conversion from magnitude to intensity). 

Our analysis adopted the interpolation method 
of calculating the c ross-correlation function (CCF, 
IWhite fc Peterson! 1 19941 ). The light curves were shifted 
backwards and forwards up to 30 days against each other 
in steps of 1 day. The results of this analysis are shown in 



figure [4] 

The CCF analysis reveals a clear anti-correlation 
between the soft and hard X-rays. On the other hand, 
the optical emission is positively correlated with the soft 
X-rays. However, we also note that the optical emission 
seems to lag the soft X-rays by ~ 1-2 days, which is much 
too long to be the result of X-ray reprocessing/irradiation 
(time-scale of ~ minutes) . The optical emission, like the soft 
X-rays, is anti-correlated with hard X-rays, and the lag is 
slightly longer. However, we caution that the lag appearing 
in the CCF may not be real, but is instead an artefact 
related to the gaps in the optical data. A more regular and 
longer term optical campaign will be needed to confirm this. 



3.3 X-ray spectral analysis 

In an attempt to understand how the X-ray spectral proper- 
ties of 4U 1636-536 evolve through the different phases of the 
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Figure 3. The long-term (2005-2009) time-dependence of power 
spectra of 4U1636-536 from RXTE/ASM (upper) and Swift/BAT 
(lower) between 2005 and 2009 showing how the qu asi-periodic 
variat ion evolves. The ~46d modulation reported by iShih et al.l 
1120051 ) is marked by the dashed vertical line. 



Figure 2. Lomb-Scargle periodograms of soft and hard X-rays 
(upper), and optical B and V band data (lower) for 4U1636-536 
in 2008. The dash-dot lines indicate the 99.9% confidence level 
(see text). 

long-term variability, a broad-band spectroscopic analysis 
was performed. We used archival RXTE/PCA and HEXTE 
data bec ause of their wide spectral range between 2 and 
250 keV (jjahoda et alj|200d : iRothschild et"al]|l998h . Two 
observations, which represent two distinct epochs in the 
RXTE/ASM light-curve (at maximum and minimum) and 
sufficiently long to provide good signal-to-noise ratio, were 
selected (see figure |5J: 

• ObsID 70036-01-02-00 (MJD 54270), at the peak of an 
oscillation 

• ObsID 70034-01-01-000 (MJD 53612), at the minimum 
of an oscillation 

The spectra were derived from PCA Standard-2 and 
HEXTE Archive Mode following standard procedures: using 
only PCU 2 (the only PCU that operates in all our datasets) 
in PCA and Cluster B (which provides the necessary back- 
ground infbjmatirjn^ 

data (|Shaposhnikov. Jahoda. &s Markwardtl 120091 1 . The 
combined PCA and HEXTE data were then analysed 
with Xspec v. 12.x. Throughout the analysis, the Galactic 
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Figure 4. Cross-correlation functions of the three 4U1636-536 
datasets: optical to soft X-rays (dot-dashed line); optical to hard 
X-rays (dashed line); and soft to hard X-rays (solid line). 
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Figure 5. (Upper) The complete RXTE/ASM light-curves of 
4U1636-536 covering more than 13 years. The data have been 
rebinned with an interval of 5 days for clarity at this scaling. 
The black triangles indicate times of RXTE observations selected 
for detailed spectral analysis (see text). (Below) Details from the 
original one day average RXTE /ASM data which are in the vicin- 
ity of the RXTE/PCA observations. 



colum n density (Nh ) toward the source was fixed at 3xl0 21 
cm " 2 (jDickev fc LockmarJll990h . We note that this value is 
well within the range determined by the ROSAT P osition 
Sensi tive Proportional Counter observations (PSPC. ISchulj 
1999). To estimate the luminosity and radii (see context 
below), the distance of 4U 1636-535 is referred to ~ 5.9 kpc 
jCornelisse et al. 1120031 ). 

Figure [5] shows the data and model fits, with table 
[1] summarising the best fit model details. The spectral 
properties of 70036-01-02-00 (RXTE/ASM high epoch) 
can be simply described by a model consisting of a multi- 
temperature accretion disc (DiskBB) component plus a 
thermal Comptonisation (CompTT, a spherical geometry was 
assumed) component. 

On the other hand, the spectral properties of 70034-01- 
01-000 (RXTE/ASM low epoch) is more complicated. The 
same two-component model fit to the data is poor (xt 3> 
2). Instead, it is more compatible with a single-temperature 
blackbody component (Bbody) plus a thermal Comptonised 
one (using the same input soft photon temperature To as 
in the high-state epoch). The blackbody component likely 
represents the emission from the surface of the neutron star. 
Additionally, a power-law component for the hard X-rays 
> 50 keV and a Lorentz profile at ~ 6.3 keV (which may 
represent a relativistically broadened Fe Ka emission line) 
are required to improve the fitting result. We note that these 
additional components are quite similar to those observed 
in the hard st ate of black hole X-ray binaries , such as GRO 
J1655-40 (see iRemillard fe McClintockll2006t ). 



ObsID 


70036-01-02-00 


70034-01-01-000 


ASM Epoch 


high 


low 


DiskBB 


T in (keV) 
nor. 


0.65±0.06 
296.5±i| 9 3 3 




Blackbody 


kT (keV) 




1.14±0.25 


nor. 




0.002±0.001 


compTT 


To (keV) 
kT (keV) 

T 

nor 


1.05±0.08 
2.88±§;o| 

c lfi+0.28 
°' io=t 0.34 

0.05±0.01 


1.05 (fixed) 

tr iciO.90 
°- io=c 0.85 

>5.13 
n nn7+ - 002 

U.UU<± 004 


Power-Law 


r 




2.12±0.04 


nor. 




0.29±°;° 3 


Lorentz 


LineE (keV) 
Width (keV) 
nor. 




6.28±°;ff 
3.74±i;i 
0.008±°;03 7 


L X (10 36 erg s" 1 ) 


2.0-80.0 keV 


4.41±0.02 


7.50±0.01 


xl 


1.20 


0.98 



Table 1. Spectral model fits to two distinct states of 4U1636- 
536 with the RXTE PCA/HEXTE data. The X-ray luminosity is 
determined from the unabsorbed flux. 



4 DISCUSSION 
4.1 X-ray variability 

From the 3 years of simultaneous RXTE/ASM and 
Swift/BAT archival data on 4U1636-536, it is clear that 
the variability, on both short and long time-scales, is 
indeed strongly energy-dependent. This was first noted 
in simultaneous R XTE/ASM and INTEGRAL/IBIS data 
ijShih et al.l l2005h . and has been a persistent feature of 
this source since 2002, when RXTE/ASM observed the 
beginning of a dramatic decline in flux. Additionally, we 
find that the long-term variability of 4U 1636-536 has been 
constantly evolving over the last decade, although not in a 
periodic manner. 

A tilted/warped, precessing accretion disc was sug- 
gested as an ex planation for such long term variations 
|Shih et al.ll2005l ), and the system parameters (in partic- 
ular the low mass-ratio) of 4U 1636-536 do allow for a 
precessing accretion disc jCasares et al.l [2006). However, 
the presence of an evolving trend in the long-term period, 
the anti-correlation between soft and hard X-rays, and the 
distinct X-ray spectral properties at different epochs are all 
difficult to understand in terms of a precessing disc. 

The X-ray spectroscopic analysis clearly reveals two 
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Figure 6. Spectra of a bright (left) and low (right) epoch observed by RXTE/PCA and HEXTE, shown together with the total model 
and residuals with respect to the corresponding best fits. 



different spectral states in the high and low epochs. It 
suggests that, in the high epoch, 4U1636-536 is like a typical 
LMXB: the accretion disc is fully developed and extends 
into the inner boundary region (Ri„ ~ 25 km, assuming 
face-on orientation, i=0°), whic h is likely surrounded by a 
Comptonised region, or corona (|Mitsuda et al.lll989l ). The 
energy state is softer as the spectrum is dominated by soft 
X-rays (2-10 keV). 

In the low epoch, on the other hand, no reliable multi- 
temperature disc component could be derived suggesting 
that the accretion disc probably retreated away from the 
inner region, or was truncated. Instead, we have determined 
a soft X-ray emitting region of radius ~ 5 km, which most 
likely corresponds to the neutron stellar surface. We also 
note that the spectral data in the low epoch contains a 
high energy component (>50 keV) which is very similar to 
that in a "peculiar" INTEGRAL observation of 4U1636-536 
l|Fiocchi et al.l [2006), as well as to those of black hole 
binaries in hard X-ray states. Clearly, 4U 1636-536 became 
significantly harder during the low epoch, as the hard X-ray 
luminosity (L x :20-sofceV ~ 3.2 x 10 36 erg s _1 ) composed 
roughly half of the bolometric X-ray flux (2.0-80 keV, see 
table [p. 



4.2 Comparison with X-ray transients 

The observed X-ray variability phenomena from 4U1636- 
536 resemble the main X-ray spectral states of black hole 
X-ray binaries, namely the "hard" and "thermal" state s 
i|Done fc Gierlinskil 120031 ; iRemillard fc McClintockl l2006h . 
Currently, various flavours of accretion flow models have 
been developed, mainly for black hole X-ray binaries, to not 
only explain these main spectral states, but al so the more 
subtl e states observed (e.g. see the review by iDone et al.l 
120071 ). The most commonly used explanation assumes that 
the optically thick inner disc is replaced by a hot, optically 
thin, flow when black hole X-ray binaries change to the 
"hard" state. Although a similar scenario is used to explain 



the state transition in atoll sources, there is a fundamental 
difference with their black hole cousins: the presence of 
an observable solid surface at the neutron star. With its 
frequent state changes and abundant archival data from 
RXTE, Swift and INTEGRAL, 4U1636-536 could be an 
ideal target to study the behaviour of the accretion flow in 
neutron star binaries, and thereby adapt the models that 
are currently used for the BHXRBs. 



4.3 Optical/X-ray correlation 

Our three-month optical monitoring campaign on 4U1636- 
536 reveals a clear modulation which was obviously 
correlated with the X-ray states. Not surprisingly, this 
is consistent with the currently accepted mechanism for 
LMXB optical counterparts that they are dominated 
by X-ray reprocessing on the accretion disc and heated 
face of the donor star by t he cen tr al X-ray source (e.g. 
van Paradiis fc McClintockl Il994l ; IShahbaz fc Kuulkersl 
1998h . Thus the observed optical component largely follows 
variations in the intrinsic X-ray flux. 

Our campaign also reveals that the increasing hard 
X-rays (> 10 keV) does not seem to influence the optical 
modulation. This is in contrast to LMC X-2 where its 
optical/X-ray correlatio n was much stronger w ith the hard 
X-rays than the soft (|McGowan et al.l 120031 ). suggesting 
that different accretion mechanisms are working in the 
inner regions of these systems. 

One possible scenario that could explain our ob- 
servation of 4U 1636-536 is Compton reflection (see 
Magdziarz fc Zdziarskilll995l ; iRevnolds fc Nowakll2003l ). In 
this case, the softer part of X-ray continuum (power-law 
component with T ~ 2.12) are photoelectrically absorbed 
by material in the disc/stellar surface, while a large fraction 
of the hard X-rays are scattered out of the system with little 
change in photon energy. Since this would not explain the 
observations of LMC X-2, we think that further research into 
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this topic is warranted, especially given the improvement in 
hard X-ray spectroscopic sensitivity over the last decade. 



5 CONCLUSIONS 

We have demonstrated that 4U1636-536 exhibits a coherent 
long-term variation across a wide range of wavelengths. 
The soft and hard X-rays are clearly anti-correlated on 
the ~30-50d variability timescale, and this appears to be 
a persistent phenomenon. The evolution of the broad-band 
spectral characteristics of 4U 1636-536 resembles the hard 
and soft/thermal state change usually seen in black hole 
X-ray binaries. 

Our work suggests that 4U 1636-536 is a unique and 
valuable target to understand the physics of X-rays binaries, 
including the size and location of the hard component in 
neutron star LMXBs, the interaction of hard X-rays with 
their surrounding environments such as the accretion 
disc/companion star surface, and the similarity or dissim- 
ilarity of accretion disc evolution between neutron star 
and black hole X-ray binaries. We are planning to explore 
these topics through a comprehensive optical spectroscopic 
and broad-band X-ray spectral analysis covering different 
variation epochs in the future. 



6 ACKNOWLEDGMENTS 

We acknowledge the quick-look results provided by the 
ASM/RXTE teanfl, and the Swift/BAT transient monitor 
results provided by the Swift/BAT tearrQ. ICS would like to 
thank the Department of Physics and Astronomy, Michigan 
State University for providing the support and accessing 
the 1.3m telescope managed by SMARTS Consortium. 
RC acknowledges a Ramon y Cajal fellowship (RYC-2007- 
01046) and a Marie Curie European Reintegration Grant 
(PERG04-GA-2008-239142). We would also like to thank 
for the anonymous reviewer's valuable comments. 



Cornelisse, R., et al. 2003, A&A, 405, 1033 
Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 
Done, C, Gierlihski, M., & Kubota, A. 2007, AAPR, 15, 1 
Done, C, & Gierlihski, M. 2003, MNRAS, 342, 1041 
Farrell, S. A., Barret, D., & Skinner, G. K. 2009, MNRAS, 
393, 139 

Fiocchi, M., Bazzano, A., Ubertini, P., & Jean, P. 2006, 

ApJ, 651, 416 
Hasinger, G. & van der Klis, M. 1989, A&A, 225, 79 
Jahoda, K., Markwardt, C. B., Radeva, Y., Rots, A. H., 

Stark, M. J., Swank, J. H., Strohmayer, T. E., & Zhang, 

W. 2006, ApJS, 163, 401 & Yamaoka, K. 1998, PASJ, 50, 

667 

Levine, A.M., Bradt, H., Cui, W., Jernigan, J.G., Morgan, 
E.H., Remillard R., Shirey, R.E. & Smith, D.A. 1996, ApJ, 
469, 33 

Liu, Q. Z., van Paradijs, J., & van den Heuvel, E. P. J. 

2007, AAP, 469, 807 
Magdziarz, P., & Zdziarski, A. A. 1995, MNRAS, 273, 837 
McGowan, K. E., Charles, P. A., O'Donoghue, D., & Smale, 

A. P. 2003, MNRAS, 345, 1039 
Mitsuda, K., Inoue, H., Nakamura, N., & Tanaka, Y. 1989, 

PASJ, 41, 97 

Reynolds, C. S., & Nowak, M. A. 2003, Physics Reports, 
377, 389 

Remillard, R. A., & McClintock, J. E. 2006, ARA&A, 44, 
49 

Rothschild, R. E., et al. 1998, ApJ, 496, 538 

Schulz, N. S. 1999, ApJ, 511, 304 

Shahbaz, T., & Kuulkers, E. 1998, MNRAAS, 295, LI 

Shaposhni kov, N., Jahoda, K., fc M arkwar dt, C. B., 2009, 

http://www.universe.nasa.gov/xrays/programs/rxte/pca/ 

doc /rmf / pcarmf- 11.7/ 
Shih, I. C, Bird, A. J., Charles, P. A., Cornelisse, R., & 

Tiramani, D. 2005, MNRAS, 361, 602 
Smale, A.P. & Mukai, K., 1988, MNRAS, 231, 663 
van Paradijs, J., & McClintock, J. E. 1994, AAP, 290, 133 
White, R. J., & Peterson, B. M. 1994, PASP, 106, 879 



REFERENCES 

Barthelmy, S. D., et al. 2005, Space Science Reviews, 120, 
143 

Belloni, T., Homan, J., Motta, S., Ratti, E., & Mendez, M. 
2007, MNRAS, 379, 247 

Casares, J., Cornelisse, R., Steeghs, D., Charles, P. A., 
Hynes, R. I., O'Brien, K. & Strohmayer, T. E. 2006, MN- 
RAS, 373, 1235 

Charles, P., Clarkson, W., Cornelisse, R. & Shih, I.-C. 2008, 
New Astron.Rev. 51, 768 

Clarkson, W. I., Charles, P. A., Coe, M. J., Laycock, S., 
Tout, M. D., & Wilson, C. A. 2003, MNRAS, 339, 447 



4 including all those working on the ASM at MIT and at the 
Goddard Space Flight Center SOF and GOF. 

5 including all those working on the BAT at the Goddard Space 
Flight Centre and Los Alamos National Laboratory 



